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t^ ' Abstract. Self-annihilating dark matter gravitationally captured by the Sun could yield 

observable neutrino signals at current and next generation neutrino detectors. By exploiting 
such signals, neutrino detectors can probe the spin-dependent scattering of weakly interacting 
massive particles (WIMPs) with nucleons in the Sun. We describe a method how to convert 
constraints on neutrino fluxes to a limit on the WIMP-nucleon scattering cross section. In 
this method all neutrino flavors can be treated in a very similar way. We study the sensitivity 
C^ '_ of neutrino telescopes for Solar WIMP signals using vertex contained events and find that this 

detection channel is of particular importance in the search for low mass WIMPs. We obtain 
highly competitive sensitivities with all neutrino flavor channels for a Megaton sized detector 
through the application of basic spectral selection criteria. Best results are obtained with 
the electron neutrino channel. We discuss associated uncertainties and provide a procedure 
how to treat them for analyses in a consistent way. 

Keywords: Dark Matter, Solar WIMPs, WIMP-nucleon scattering cross section. Neutrino 
Astronomy 



1 Introduction 

There is compelling observational evidence for the existence of dark matter; however, its 
nature remains unknown. A prominent candidate that naturally arises in supersymmetric [1] 
and large extra dimension models [2] is a weakly interacting massive particle (WIMP) [3]. 
WIMPs are predicted to have masses ranging from a few GeV to several TeV [4]. They 
could be detected through a variety of possible signals, indirectly through the byproducts 
of self-annihilations or decays, directly by scattering effects with nuclei, or production at 
colliders. 

The nucleon coupling of a WIMP in the extreme nonrelativistic limit is characterized 
by a scalar (spin-independent) and an axial vector (spin-dependent) term. In general the 
spin-dependent (SD) couplings for WIMPs on protons and neutrons are expected to be of the 
same order of magnitude, and very similar for the spin-independent (SI) case [5]. Throughout 
the paper we assume that the WIMP-proton and WIMP-neutron scattering cross sections 
are identical and do not consider isospin violating interactions, which have been discussed 
elsewhere (see for example ref. [6]). 

The spin-dependent and spin-independent terms behave very differently if coherence 
across the nucleus is taken into account; the nucleon contributions interfere constructively 
to enhance the WIMP-nucleus elastic cross section. The WIMP-nucleon cross sections for 
spin-independent interactions are generally orders of magnitude smaller than those allowed 
for the spin-dependent case. Constructive interference across the nucleus leads to a quadratic 
enhancement with the nucleon number (~ A^) in the spin- independent couplings. By exploit- 
ing this dependence through the use of high Z elements, direct detection experiments such 
as CDMS-II [7], XENONIOO [8], and EDELWEISS-II [9] have been able to tightly constrain 
the spin-independent WIMP-nucleon scattering cross section. 

Limits on the spin-independent cross section extend to 10"'^ pb (10~^^ cm^) for WIMPs 
mass around 50 GeV. For the spin-dependent case constraints from direct detection experi- 
ments, such as SIMPLE [10], COUPP [11], PICASSO [12] are much weaker and range down 
to a minimum of about 2.8 x 10"^ pb for a WIMP mass of 45 GeV. Better sensitivities are 
achieved by indirect searches with neutrino telescopes for this case. Depending on the as- 
sumed annihilation channel, limits for the spin-dependent case from Super-Kamiokande [13] 
and IceCube [14] range from 10~^ pb to 10~^ pb. These tight constraints are obtained by 
looking for neutrino signals from self-annihilating dark matter gravitationally captured by the 
Sun. The comparison becomes possible, as the capture process is initiated by the same inter- 
action that direct detection experiments exploit: a nuclear recoil in which a WIMP transfers 
a fraction of its kinetic energy to a nucleon. For the commonly assumed case of equilibrium 
between capture rate, Fc, and annihilation rate. Fa, neutrino fluxes, cpl, from the Sun are 
directly dependent on the WIMP-nucleon scattering cross section. Hence, the neutrino flux 
from the Sun or any constraint on it can directly be linked to the WIMP-nucleon scattering 
cross section. 

While previous work has focused on the problem of converting the WIMP-nucleon scat- 
tering cross sections to muon neutrino generated muon flux, (pj^ [15, 16], we here focus on 
the neutrino flux instead. Muon fluxes inherently contain a dependence on the detector lo- 
cation as well as assumptions on neutrino interactions and muon propagation through some 
medium. Using the neutrino flux itself has the benefit that it allows to treat all neutrino 
flavors in a very similar way and the muon neutrino induced muon flux is just one special case 
of it. The neutrino flux is directly related to the study of vertex contained events at neutrino 



detectors. This class of events is in particular relevant for precision measurements of neutrino 
flavor ratios and the extraction of neutrino spectral information. We discuss results in light 
of a generic neutrino detector. Results are applicable for example to Super-Kamiokande, 
IceCube/DeepCore, KAMLAND, and future detectors proposed for long baseline neutrino 
oscillation experiments such as Hyper-K [17] and LENA [18]. For a detector size (diameter) 
that is on the same order or larger compared to the average length of a muon track created 
by a muon neutrino with energy of interest, vertex contained events will be the dominant 
part of the expected signals in the detector. 

This paper is structured as follows: In section 2 we provide a description how to compare 
a neutrino flux from the Sun to the spin-dependent WIMP-nucleon scattering cross section, 
a . Our result is based on the DarkSUSY package version 5.0.4 [19], which allows one to 
obtain conversion factors for given halo models and supersymmetric model parameters. Our 
study is focused on the WIMP mass range favored by supersymmetry (10 GeV to 10 TeV), 
however results are completely model-independent and applicable to any model that produces 
WIMPs in this mass range. We discuss the prospects of observing WIMP signals at neutrino 
detectors in section 3. We use a general approach that utilizes vertex contained events of all 
active neutrino flavors. Results are quoted for a data size of 200 kton-years and 5 Mton-years. 
We calculate differential neutrino fluxes, compared to detector backgrounds from atmospheric 
neutrinos and give event spectra. In section 4 we describe how to incorporate and treat 
uncertainties related to the flux conversion. In particular we address uncertainties due to 
the dark matter distribution, which is underlying to all searches. We study the equilibrium 
condition and discuss procedures for scenarios when it is not satisfied. We summarize and 
draw our conclusions in section 5. 

2 Neutrino Flux Conversion 

In this section we review the neutrino flux conversion to the WIMP-nucleon scattering cross 
section as implemented in DarkSUSY [19]. We introduce the parameters and describe the 
choice of values used in this study. 

2.1 Revievif of Standard Calculation 

WIMPs from the Galactic dark matter halo could be gravitationally captured by the Sun. 
This process is initiated by an initial scatter of a WIMP with a nucleon in the Sun, in which 
the WIMP could lose enough energy to be gravitationally bound and eventually captured 
by the Sun. The number of WIMPs in the Sun, M, is governed by the following differential 
equation with three constants, 

^ = Cc - CaAA2 - CeAT. (2.1) 

at 

Cc is a WIMP capture term, that depends on the WIMP-nucleon scattering cross section. 
It can be identified directly with the capture rate (Cc = Pc). Ca controls annihilation and 
depends on the dark matter self-annihilation cross section, (aAv). Ce is an evaporation term, 
which is only relevant for WIMP masses below 4 GeV [20, 21], and is therefore outside the 
mass range considered here. Neglecting Ce, the annihilation rate. Fa, is given by Fa = 
i^Cp^M"^. As the Sun accumulates dark matter the annihilation rate at time t can be linked 
to the capture rate via the annihilation equation. Under the assumption of a static capture 
rate 



rA = ^Cctanh2(t/r), (2.2) 



with the equihbrium time scale, r, given by 



/(-)/- = A/T^r^{exp(-^^^^— f^ -expp:--^)), (2.5) 



r = 1/VCcCl. (2.3) 

For the case that t ^ t equihbrium between capture and annihilation is reached. At present, 
using the age of the Sun with t = t0 ~ 4.7 x 10^ years the equilibrium condition is in general 
fulfilled.^ For this case the annihilation rate is regulated by the capture rate, Fc and hence 
independent from the dark matter self-annihilation cross section: 

Ta = Icc. (2.4) 

We follow earlier approaches [16, 22] as implemented in DarkSUSY for the calculation of the 
capture rate. For the velocity of the WIMPs outside the potential well of the Sun in the 
Galactic frame, u, we assume a Maxwellian distribution, 

3iu-VQf \ _ ^^ ( SJu + vq ) 

where Vq is the velocity of the Sun relative to the halo, Vd is the WIMP velocity dispersion, 
and n^ is the local WIMP number density given by the ratio of local density, po, and WIMP 
mass, ruy.. Using this assumption on the velocity, the total capture rate can be written as, 

Cc = lJ"^^r^drY^^, (2.6) 

where Rq is the radius of the Sun, Ci is the WIMP capture rate per unit shell volume of 
given composition from element "z" in the Sun. We use the solar composition model BS2005- 
AGS [23]. Since the Sun is predominately a Hydrogen target WIMP-nucleon scatterings are 
expected to be dominated by spin-dependent interactions. The capture total rate Fc is the 
sum of a spin-independent F^ and spin-dependent F^^ term. We assume F^ 3> Fq and 
neglect spin-independent interactions in the Sun by setting the spin-independent WIMP- 
proton cross section to zero, so that Fc ~ ^c^- This choice is well motivated by the very 
stringent constrains on a from direct detection experiments [7-9] . Alternatively, one can use 
a neutrino flux limit from the Sun to probe the spin-independent cross section by assuming 
Fc ~ Fq . For this condition in practice results are roughly two orders of magnitude smaller 
compared to the spin-dependent scattering cross section (see e.g. ref. [24]). 

The WIMP capture rate per unit shell volume from element i in the Sun is given [22] 
by 

dCi r— f{u) 

—— = I du wil^,i{w), (2.7) 

dV Jo u ' 

where at the interaction point in the Sun, the WIMP has been accelerated by the corre- 
sponding escape velocity, v, so that its velocity is given hy w = ^/v? + v"^ . Qi,^i{w) is the 
capture probability per unit time for the corresponding element i (we refer the reader to 
Wikstrom and Edsjo [16]). The upper integration bound Umax is the velocity at which the 



^We discuss a treatment for deviations from the equilibrium between capture and annihilation in section 4.1. 



Table 1. Capture rate as function of WIMP mass and scattering cross section a^^ as obtained with 
DarkSUSY and compared with resuhs from the analytic approximation eq. (2.8) [25, 26]. We show the 
impact of inchiding and neglecting planetary effects following the default DarkSUSY implementation 
of them. 



WIMP mass 


Capture Rate Tq / cr 


^^ [(l/s)/pb] 


^x 


Jupiter 


Approximation 


(GeV) 


neglected 


included 


eq. (2.8) [25, 26] 


10 


9.43 X 10^'^ 


9.42 X 10^'^ 


2.25 X 10^^ 


50 


8.74 X 10^6 


8.69 X 10^6 


9.00 X 10^6 


100 


2.48 X 10^6 


2.45 X 10^6 


2.25 X 10^6 


500 


1.09 X 10^5 


1.01 X 10^5 


0.90 X 10^5 


1000 


2.75 X 10^4 


2.38 X 10^^ 


2.25 X 10^4 


5000 


1.11 X 10^3 


4.82 X 10^2 


9.00 X 10^2 


10000 


2.77 X 10^2 


3.47 X 10^1 


2.25 X 10^2 



WIMP scatter to the escape velocity [16]. The rough dependence of the capture rate on 
various quantities can be understood from the following analytical expression [25, 26]: 



24„-l 



Cc ^ 9 X lO^'-'s 
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10-2pb 



50 GeV 



m. 



(2., 



While equation (2.8) is good to demonstrate the dependencies, rates deviate significantly 
compared to the full computation with DarkSUSY, as shown in table 1, and should not be 
used for event rate estimates. 

The differential neutrino {vj + Uj) energy spectrum of flavor j at Earth per annihilation 
is given by: 



n 



f 






E 



f 



dm 



f 



dE^^ 



(2.9) 



where P(j, i) is the probability that produced neutrino of flavor i oscillates to flavor j when 
they propagate to the detector, Bf is the branching ratio for annihilation channel /, and 
dNJ^ /dEy. the corresponding differential number of neutrinos of flavor i per annihilation. 



The integrated neutrino flux for an annihilation rate. Fa, above an energy threshold E is 
then given by 



^L 



dN, 



E>E^ 



dAdt 



Ta 



^th, 



dEy^nl^. 



(2.10) 



Here Dq is the distance of the center of the Sun to Earth. 

For our study we consider several benchmark annihilation channels that represent soft 
(66) and hard (VF~'~VF~,r''~r~) neutrino spectra. We also take into account the suppression 
of the capture rate due to the gravitational effects from planets, especially Jupiter, during 
the capture mechanism [27, 28]. We use the default DarkSUSY implementation of this 
process. We scan through MSSM-7 models [29] using DarkSUSY and enforce the equilibrium 
condition, so that the flux is independent of the evolution of the dark matter in the Sun. 
We calculate conversion factors and the differential neutrino flux from WIMP annihilations 
based on the assumptions and equations above. 



Oh 



'b 10 




10' 10" 

m^ (GeV) 



10" 



Oh 



10' 



10' 



10' 



10' 



10' 



D 10- 



10' 



■11 



■12 



■13 



■14 



■15 



■16 



■17 



10" 



' KRv '....^.... 


1 i 


DD V^ 


/ 


tt V- 


y^-'A 




yy ■ 


T X V^ 


/r ' 


' y 


y ■ 


" y^ 


1 


r .••''' / 


n 


f"' ,....-■'■'" 


1 




1 " 



10' 



10' 



10^ 



10" 



nij^ (GeV) 



Figure 1. Conversion factor, which we computed with DarkSUSY, for the integrated neutrino flux 
at Earth above a threshold of 1 GeV to the spin-dependent WIMP-nucleon scattering cross sec- 
tion obtained through DarkSUSY. Left: Neutrino flux conversion for different neutrino flavors to 
^SD £qj, g^ given annihilation channel (66). Right: Flux conversion for various annihilation channels 
(66, ti, VF"'"VF~,r"'"T~) for tau neutrino fluxes. 



At production the neutrino flavor ratio is given by [30]: 



Vp-.Vp-. v., 



1 : 1 : 1 : 1 :r 



(2.111 



where the factor, r, depends on the annihilation channel. Oscillation and propagation effects 
distort this ratio, so that at Earth the yields for muon and tau neutrinos are more or less equal. 
All oscillation and propagation effects were taken into account as part of DarkSUSY [19]. 
For the neutrino oscillations, we used the parameters described in ref. [31] as dyi = 33.2° , 
6li3 = 0°, 6*23 = 45°, 5 = 0, Am|i = S.lxlO-^eV^ and |Am|i| = 2.2'Kl{)-^e\'^ . For aneutrino 
energy threshold of E^"^ = 1 GeV, we compute conversion factors for all three neutrino 
flavors and various annihilation channels. We assume that the local WIMP density pQ is 
0.3 GeV/cm'^, a Maxwellian velocity distribution with a velocity dispersion of f^ = 270 km/s. 
Figure 1 illustrates conversion factors for various neutrino flavors and annihilation chan- 
nels that we computed with DarkSUSY. a^^ /(pi is nearly identical for z/^ and u-j- and the 
difference relative to u^ is roughly 20%-30% driven by differences in the injection spectrum 
and oscillations. 



3 Results 

In this section we discuss the differential neutrino flux from WIMP annihilations in the 
Sun and compute event rates at detectors. We assume a generic detector and quote results 
per kiloton fiducial volume. We do not consider events interacting outside the detector 
volume, which have been studied extensively in previous works [13-16], and focus instead 
on vertex contained events. The usage of vertex contained events for analyses has been 
well established and is documented for Super-Kamiokande [32]. Rates for these events are 
easily scalable to any detector size and results are to first order independent of the actual 



detector geometry. Rates are given as differential rates as function of energy and would have 
to be convoluted with the energy response and resolution for a specific detector. We chose 
a coarse energy binning to mimic a limited energy resolution. To put results in context, we 
compare event rates with those from atmospheric neutrinos [33]. We use the angle-averaged 
flux and neglect atmospheric neutrino oscillations as these are most relevant at an energy 
range, which is largely below our region of interest and in addition results would strongly 
dependent on the detector location. The atmospheric tau neutrino flux component, which 
is almost entirely generated through neutrino oscillations, has a very strong dependence on 
the underlying baseline length and therefore zenith angle. We compare signal muon and 
electron neutrino rates to the corresponding component from the atmospheric neutrinos. As 
tau events are in general more similar in topology to electron neutrino events, we compare 
the tau neutrino signal rates to those from atmospheric electron neutrinos. We assume the 
atmospheric neutrino background is well determined, which is a reasonable assumption as it 
can be estimated from the data itself from off-source regions. Further the size of the available 
off-source region is large compared to the on-source region, which minimizes the uncertainty 
on the background prediction. We compare the sensitivity of different search strategies for a 
generic neutrino detector assuming an angle-averaged atmospheric neutrino flux. We discuss 
other Sun related backgrounds in section 3.3. 

3.1 Expected Event Rates 

We compute event rates for charged-current (CC) interactions, as they allow flavor iden- 
tification and measurement of the neutrino energy. We do not consider the impact of 
neutral-current (NC) interactions in this work. Deep inelastic neutrino and anti-neutrino 
cross sections are treated separately [34] and we use the ratios of o-'^%j/o'^%s to treat tau 
neutrinos [35, 36]. 

We order results by neutrino fiavor and compare the rates to those produced by the 
corresponding angular averaged atmospheric neutrino component. The differential event 
rate at the detector for each channel is be given by 

'-^ = T--^<''N{E)niiE)e., (3.1) 

where n is the number of nucleons per kton given by Na x kton ~ 6.023 x 10^^, where Na is 
the Avogadro number. T is the detector livetime in units of kton-years. nl^ is the neutrino 
spectrum at Earth per annihilation computed with DarkSUSY as defined by equation (2.9). 
As annihilations occur almost entirely near the center of the Sun, the resulting neutrino flux 
can be considered a point-like source. The acceptance, e^, deflnes the fraction of events that 
will be contained within an opening angle around the Sun. 

We assume three different options for the opening angle around the Sun and do not 
optimize them to keep results from being exposure-dependent. We use two fixed angle cuts 
of 30° and 10° and a neutrino energy dependent selection criteria ^{E), which is used as 
our benchmark cut. The opening angle contains 68% of the signal under the assumption 
of the kinematic angle between neutrino and lepton convoluted with an intrinsic detector 
resolution of one degree. Hence, ^(i^jy) — l°y^l + YTeV/E^. This case requires that the 
energy of vertex contained events can be reasonably well measured on an event by event 
basis to extract spectral information. 

In figure 2 we compare the differential neutrino event rates to those caused by atmo- 
spheric neutrinos. We use a WIMP-nucleon scattering cross section of 1.0 x 10~^ pb, which is 



of the same order as present exclusion limits [13, 14]. For the corresponding annihilation rate 
please refer to table 1. We show the rates as function of neutrino energy for neutrinos and 
anti-neutrinos separately and compare results with the atmospheric neutrino background. 
The number of background events depends on the angular opening angle around the Sun; 
hence the background is given by: 



dE 



Tna%l{E)^tl^{E) 



^(E) 



2-Kcos{e)de. 



(3.2) 



For cascade-like events the rate would therefore correspond roughly to the number of 
fully contained events, while for the muon neutrino flux the given rate corresponds to the 
fully and partially contained events. We have cross-checked our estimates against reported 
multi-GeV event rates of the corresponding types at Super-Kamiokande [37] and found good 
agreements. 
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Figure 2. Rates for neutrino interactions (z/o + i'c), v^-, and i7^ per kton-years for a solar WIMP 
signal for various annihilation channels compared to those from atmospheric neutrinos. We assume 
an opening angle selection criteria containing 68% of the signal. The number of events per year is 
given for each energy bin. The neutrino channel for tti^ = 10 GeV (left) and ttI;^ = 1 TeV (right) are 
shown for v^.. A spin-dependent WIMP-nucleon scattering cross section of 1.0 x 10~^ pb (1 fb) was 
assumed, which leads to an annihilation rate of 9.42 x lO^^'s"^, and 2.38 x lO^-'^s^^ at equilibrium for 
Wimp masses of 10 GeV and 1 TcV, respectively. 

Since neutrino detectors are generally not sensitive to charge, except magnetized designs 
like MINOS [38] , we give rates for the sum of neutrinos and anti-neutrinos. Figure 3 shows 
event rates per kton-year for CC interactions by neutrino flavor compared to the correspond- 
ing background, which would result in track or cascade-like events. As an example we have 
calculated event rates for 200 kton-years, which roughly corresponds to the total integrated 
lifetime of 14 years of Super-Kamiokande data, see table 2 for event rates assuming a lower 
energy threshold of 1 GeV. Hyper-K or similar next generation water Cherenkov detectors 
would acquire a comparable lifetime in a year or less [17, 39]. Rates need to be compared 
with about 300 (646) atmospheric electron (muon) neutrino events for a threshold energy of 
1 GeV. For energies above 10 GeV rates drop to 2 (8) for electron (muon) neutrinos. 

3.2 Sensitivity 

We compute the sensitivity for a generic neutrino detector assuming a data sample, of 
200 kton-years and 5 Mton-years, as could be achieved by a future detector such as Hyper- 
K [17], LENA [18], or an infill array to IceCube/DeepCore. For the sensitivity study we com- 
pare the expected neutrino spectrum from background only with that containing in addition 
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Figure 3. Rates per kton-years for neutrino interactions (v + D) from atmospheric neutrinos are 
compared to those from a solar WIMP signal for various annihilation channels. Neutrinos are selected 
within an opening angle around that Sun contains 68% of the WIMP signal. The event rates are 
given per energy bin for the three active neutrino flavor channels i^e, i^^, and v^ (top to bottom) and 
for two different WIMP masses m^ = 100 GcV(lcft) and m^ ~ 1 TeV (right). A spin-dependent 
WIMP-nucleon scattering cross section of 1.0 x 10^'^ pb (1 fb) was assumed, which leads to a capture 
rate of 2.45 x 10^^ 1/s (m^ = 100 GeV) and 2.38 x lO^i 1/s (m^ = 1 TeV). 



to the background a signal. We assume that the WIMP mass and dominant annihilation 
channel is unknown, and hence we consider the entire energy spectrum. In general one can 
assume that a signal would be clustered around a few adjacent energy bins. A sliding energy 
window would in this case be an optimal search strategy, but it inherently also generates 
trial factors that need to be considered. As contributions to the background are tiny from 
high energy events due to the steeply falling atmospheric neutrino spectrum, we consider an 
energy threshold rather than a window. For simplicity we use three different energy thresh- 
olds: 1 GeV for m^ < 100 GeV, 10 GeV for 1 TeV >m^> 100 GeV, and 100 GeV threshold 
for m^ > 1 TeV. As we expect significant absorption in the Sun at energies close to 1 TeV, 
moving the threshold any higher is not desired. 



Table 2. Integrated number of events in 200 kton-years. assuming a half cone opening angle containing 
68% of the signal. Using a lower energy threshold of 1 GeV we have integrated Ve,v^, v-r over an energy 
range of 1 GeV to 10 TcV. The total rate is broken down in the that from neutrinos and anti-neutrinos 
separately in the bracket. A spin-dependent WIMP-nucleon scattering cross section of 1.0 x 10^'^ pb 
was assumed, which leads to a capture rate of 2.45 x lO^'^ 1/s (m^ = 100 GeV) and 2.38 x 10^^ 1/s 
(m^ = 1 TeV). 



Flavor 
channel 



WIMP mass 
(GeV) 



Events per 200 kton-years 
XX -^bb XX ^ T+T- XX -^ W+W" 



10 

100 

1000 



9.1 (5.5/3.5) 

1.5 (0.9/0.6) 

0.05 (0.03/0.02) 



67.1 (45.8/21.3) 

12.4 (7.9/ 4.5) 

0.28 (0.15/ 0.13) 



0.0 (0.0/0.0) 

5.3 (3.3/2.0) 

0.08 (0.04/0.04) 



10 

100 

1000 



12.5 (8.8/3.7) 

1.9 (1.3/0.6) 

0.05 (0.03/0.02) 
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We determine the sensitivity based on the average 90%C.L. Feldman-Cousins [40] upper 
limit that would be obtained by an ensemble of experiments with no signal in the presence of 
a known background [41]. For large statistics (number of background events above 100), we 
approximate the expected limit with Gaussian statistics via v6x 1.645. The assumption that 
the background is well known is valid for our case as it can be determined from the data itself 
using an off-source region which is significantly larger than the target region. As the size of 
the signal depends directly on the WIMP-nucleon scattering cross section, we determine the 
corresponding a for the limit on the expected event rate. The obtained sensitivities are 
shown in figure 4. We compare three different methods: A fixed 30° cone centered on the 
Sun, a fixed 10° cone, and an adaptive cone that contains 68% of the signal for each energy 
bin as described above. Using a small or adaptive cone size generally performs best. 

Assuming a comparable angular resolution is achieved for all neutrino flavors, the elec- 
tron neutrino channel generally outperforms other flavor channels. The tau neutrino channel 
is also very compelling for WIMP masses above 100 GeV, while at lower WIMP masses rates 
are significantly decreased due to the smaller tau neutrino cross section [35, 36]. Tau events 
can nevertheless be important as they could be detected nearly background free based on 
their unique topology. We have neglected NC induced events, which will further increase 
signal rates. 

3.3 Discussions 

Besides the atmospheric neutrino background. Solar WIMP searches are also prone to irre- 
ducible neutrino backgrounds from the Sun itself. Interactions of cosmic-rays with the solar 
atmosphere have been predicted to produce 7-rays and high energy neutrinos [42-44]. These 
albedo neutrinos are produced through the decay of charged pions in the solar atmosphere. 
Flux predictions differ widely and are beyond the scope of this paper. Hence, we only give the 
magnitude of this flux, which is estimated to be about 0.01-0.04 events per year per Mton for 
each neutrino flavor assuming an energy threshold of 10 GeV [45]. These rates are however 
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Figure 4. Top left: Sensitivity for a data size of 200 kton-years for annihilations into t^t~ and 
&&. Three different analysis methods using a fixed 30° cone centered on the Sun, 10° cone, and 
an adaptive cone that contains 68% of the signal for each energy bin. Top right: Sensitivity for 
5 Mton-years for annihilations into t'^t~ , W~^W^ , and bb compared to present bounds from direct 
and indirect detection experiments. Limits from Super-Kamiokande (SuperK) [13], IceCube [14], 
COUPP [11], and PICASSO [12] are shown. Bottom: Comparison of the different neutrino flavor 
channels for 5 Mton-years for an opening angle of 30° (left) and one containing 68% of the signal 
(right). 

still small compared to the corresponding signal rates for the achieved sensitivities. Transient 
backgrounds caused by solar flare neutrinos, which can reach energies up to a few 10 GeVs 
with a peak flux at around 100 MeV for a few minutes, generally fall below the relevant 
neutrino energies discussed in this work, but might need to be considered in analyses [46]. 

Previous searches have focused on the muon neutrino induced muon flux at neutrino 
telescopes for obvious reasons: (1) the good pointing resolution of high-energy muons allows 
us to define a small search window around the Sun. It can further be used in rejecting 
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the large background caused by down-going atmospheric muons. (2) The long range of high- 
energy muons (for example ~ 1 km at 200 GeV in water /ice) , allows to detect a flux of muons 
caused by i/^ interactions outside the instrumented detector volume. Advances in detection 
methods with operating neutrino detectors and prospects for future detectors make it feasible 
to expand Solar WIMP searches to other neutrino flavors. The case of using electron neutrinos 
in Liquid Argon detectors have been discussed recently [47], and obtained a similar sensitivity 
compared to results presented here. Not only could these provide additional information 
on the WIMP properties, in case of a detected signal, but also searches relying on vertex 
contained events can generally achieve a better energy resolution. Further, if the muon range 
is of the same order as the size of the detector, then there is no longer a significant increase 
in effective area due its range. We also note that neutrinos fluxes with energies above 1 TeV 
are severely attenuated and neutrinos above 100 GeV encounter significant absorption effects 
in the Sun. Hence, neutrino searches are limited to an energy range below a few hundred 
GeV. 

4 Discussions of Sources of Uncertainty 

In this section we discuss uncertainties on the conversion of the neutrino flux to the WIMP- 
nucleon scattering cross section and provide a description for the treatment of them. 

We group the systematic uncertainties into three different categories that are uncorre- 
lated to each other: (1) uncertainty on the annihilation rate at equilibrium, r^"\ (2) un- 
certainty on the neutrino flux related to the neutrino propagation to Earth, (3) uncertainty 
related to the flux measurement at the detector and propagation through Earth. 

Through this grouping, the error treatment can be simplified as categories are inde- 
pendent of each other and can be considered separately. In addition different experimental 
results using this description can easily be compared directly. While, the first two categories 
are common to all neutrino detectors, the third category is highly detector dependent. We 
focus our discussions on the first two categories and merely provide a summary of the rele- 
vant effects for the third category. A detailed treatment would have to be performed by the 
corresponding experiments. 

For the uncertainty discussion we assume the equilibrium condition and address the 
impact of systems that are perturbed from it separately. This is part of a discussion on the 
equilibration time scale, WIMP-nucleon scattering cross section and self-annihilation cross 
section in section 4.1. 

4.1 Annihilation Rate Uncertainty 

At equilibrium the annihilation rate is solely determined by the capture rate, F^^' = gTc- 
Hence, we can use the uncertainty on Fa at equilibrium to group all uncertainties related to 
the capture and annihilation process. A variety of uncertainties need to be considered, that 
are related to the dark matter halo density and velocity distribution at the solar Galactic 
orbit, capture mechanism, and composition of the Sun. We summarize the uncertainties here 
and quantify them where possible. 

First of all, the capture rate in the Sun depends linearly on the local dark matter 
density. Hence, under the assumption of a constant capture rate as resulting from a smooth 
dark matter halo, the uncertainty on Fa is directly dependent on the uncertainty due to the 
halo density: -jr^ ~ — ^. Note that this uncertainty affects direct detection experiments 
in the same way. In our studies we have adopted a local dark matter halo density of po = 
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0.3 GeV/cm^ at the location of the Sun, which is commonly used [48]. However, recent 
measurements favor a higher value of 0.39 it 0.03 GeV/cm^ [49]. Simulations of Milky Way 
type galaxies, that include baryons, leads to a significant flattening of the dark matter halo 
in the direction normal to the stellar disk. Pato et al. [50] found that this effect could 
lead to a dark matter overdensity in the local disk of up to 41% and estimate an average 
enhancement of 21% for such a triaxial profile compared to a spherically symmetric Einasto 
profile, resulting in po = 0.466±0.033(stat)ib0.077(syst)GeV/cm^. Based on various methods 
a general be bound can be found for 0.25 GeV/cm^ < pQ < 0.70 GeV/cm^ [51]. 

The capture rate depends in a complex manner on the circular velocity of the Sun 
(relative to the halo) and the WIMP velocity dispersion, which is commonly assumed to be 
Maxwellian as given by equation (2.5). We find that the uncertainties related to the velocity 
dispersion and circular velocity is decreased for smaller WIMP masses. Figure 5 shows the 
change in capture rate for various halo parameters differing from our default Maxwellian 
velocity distribution with a dispersion of v^ = 270 km/s and velocity of the Sun relative 
to the halo is assumed to be vq = 220 km/s. Note, that recent simulations [52] indicate 
deviations from Maxwellian velocity distributions. For a discussion of the impact of different 
velocity distributions on the capture rate, we refer the reader to Bruch et al. [53], which 
have studied this effect in light of the merger history of the Milky Way. Implications for 
direct-detection experiments that are mostly sensitive to the high velocity tail have been 
discussed [54]. Contrarily for the capture rate in the Sun the low velocity region is most 
relevant. It was found that the annihilation rate in the Sun could be enhanced by as much 
as an order of magnitude relative to standard halo models if the ratio of the accreted dark 
matter is of the same order as the halo dark matter density. This increase owes to the higher 
phase space density at low velocities in the accreted dark disc. 
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Figure 5. Left: Change in capture rate as function of WIMP mass and halo parameters compared 
to the nominal model assuming a circular velocity of 220 km/s and a velocity dispersion of 270 km/s. 
Right: Impact of Jupiter on the capture rate as function of the WIMP mass. 

Galactic substructure, or subhalos, would naturally lead to varying dark matter density 
and velocity distributions that can be address with the same error treatment described above, 
assuming extremes between dense and depleted regions. The assumption that the capture 
rate is constant with time is no longer true for this case and hence the annihilation rate would 
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also be perturbed from its nominal value. Extreme cases of substructure and the impact on 
capture rates and equilibration time scale has been discussed elsewhere [55]. Contributions 
from additional dark matter populations (for example in form of a dark disk [53] ) as the Milky 
Way could have acquired during its complex merger history falls in the same category and 
is not discussed in further detail. Deviations from a smooth halo profile at the solar circle 
through substructure would introduce a time dependence on the strength of the neutrino 
flux from the Sun. N-body simulations disfavor any significant sub structure on scales that 
can currently be resolved at the solar circle [52]. There are indications for a more complex 
structure in the velocity distribution. For the case of the Sun small variations in the local dark 
matter density and velocity distribution can generally be neglected compared to the overall 
uncertainty on the average dark matter density and velocity distribution, as the expected 
annihilation rate in the Sun is better defined by the dark matter capture integrated over a 
time scale similar to the equilibrium time scale. 

The dynamics of the solar system impact the capture process and can result in a change 
in the capture rate [27]. Modeling of this effect is extremely complex and has only been 
attempted for a simplified three-body system, taking the gravitational effects of Jupiter into 
account, which is expected to be the dominant effect. Based on these studies suppression 
in the capture rate was observed. Taking additional planets into account in simulations will 
further reduce the capture rates, though Jupiter is expected to have the largest impact [56]. 
The capture rate suppression shows a strong dependence on the WIMP mass and becomes 
relevant for WIMP masses in the TeV range as shown in figure 5. To take this reduction into 
account one can define a WIMP mass dependent suppression factor a{m.y^) = (Fc/r^"™"^ )~^. 
The uncertainty on this suppression factor directly propagates to the uncertainty in Fa- Due 
to the complexity of the evaluation of planetary effect the uncertainty is not well understood. 
Conservatively one could set an upper bound on the uncertainty by comparing the capture 
rates with and without the effect from Jupiter taken into account. While this might be 
practical for WIMP masses below one TeV, where the difference is less than 15%, a more 
detailed study would be needed for higher masses. 

Uncertainties in the composition of the Sun will also alter the capture rate, as well as a 
possible spin-independent component. However, the uncertainty in the elemental abundances 
in the Sun only affects the fiux predictions on the percent level [57]. 

Table 3 summarizes the systematic uncertainties on the annihilation rate. We consider 
three different cases as examples, a low WIMP mass 10 GeV, a 100 GeV WIMP, and a high 
WIMP mass 1 TeV. Uncertainties are given with respect to the assumed default parame- 
ters: p = 0.3 GeV/cm^, and a Maxwellian velocity distribution with a velocity dispersion 
of 270 km/s and a solar circular velocity of 220 km/s. While, some uncertainties can be 
large, they generally result in an increased annihilation rate, so that constraints tend to 
be conservative. The two dominant effects, the velocity distribution and local dark matter 
density affects direct detection experiments at the same magnitude. While, the local dark 
matter density affects rates at direct detection experiments in the same way, the impact of 
the velocity distribution can have very different affects. 

While the equilibration time scale is typically short compared to the age of the Sun, in 
some scenarios the Sun might not be equilibrium. In this case the annihilation rate would 
be reduced to a fraction of the maximal rate achieved at equilibrium. The expected neutrino 
flux is reduced by the same factor. We determine the ratio of the expected flux compared to 
the maximal flux. The average flux fraction obtained via DarkSUSY is shown in flgure 6. It 
can be seen that the equilibrium condition is in general satisfied for the region of parameter 
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Table 3. This tabic suminarizes the systematic uncertainties on the annihilation rate AFa/Fa 
at equilibrium with respect to the used default values and their uncertainty. For details see text 
and references. This summary table is meant to give a rough idea on the uncertainties and their 
approximate size, an analysis would require a more detailed discussion of them including considering 
correlations between some of the quantities. 
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Figure 6. Left: Average factor by which the age of the Sun exceeds the equilibration time scale, as 
function of the self-annihilation cross section and the WIMP-nuclcon spin-dependent scattering for 
a WIMP of mass 10 GeV. Right: The corresponding fraction of the neutrino flux compared to the 
maximal flux reached at equilibrium. The dotted line indicates the thermal relic cross section, while 
the dashed lines indicate the affect on a constrained on the WIMP Nucleon cross section if the dark 
matter self-annihilation cross section is reduced so that the Sun is not in equilibrium anymore. 



For a detailed discussion on astrophysical uncertainties we would also like to refer the 
reader to Serpico and Bertone [58]. 
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4.2 Neutrino Propagation Uncertainty 

The second class of uncertainties, which can be treated independently from the uncertainty 
on the annihilation rate, Fa, is related to the propagation of neutrinos from the center of the 
Sun to Earth. 

The uncertainties on the neutrino flux, cj)i, are related to neutrino oscillation and ab- 
sorption effects as well as regeneration of (tau) neutrinos in the Sun. These effects are fully 
implemented in DarkSUSY. For a general discussion on the neutrino fluxes [30] from dark 
matter annihilation and their flavor ratios [59] we refer the reader to the associated references. 

For neutrino production at the center of the Sun, absorption is relevant for neutrinos 
above 100 GeV, and leads to a significant suppression for neutrinos above 1 TeV. The effects 
of neutrino oscillations and interactions have been studied elsewhere [60] using an event-based 
full three flavor Monte Carlo. It was concluded that for neutrinos the main effects of neutrino 
oscillations are the effective mixing of muon and tau neutrinos during the propagation to the 
solar surface and the consequent mixing of electron neutrinos during the propagation from the 
surface of the Sun to the Earth. The mixing between electron and other neutrino flavors is not 
complete, due to the leptonic mixing angle 9i2 being less than maximal, which results in the 
yields for muon and tau neutrinos to be similar but differ from those of the electron flavor. 
The impact of the oscillation and propagation effects can vary based on the annihilation 
channels and WIMP masses, and hence need to be considered for individual cases. As an 
example the systematic uncertainty due to neutrino oscillation and propagation effects was 
found to be 4% for the IceCube Solar WIMP analysis [14, 61]. Note also that the energy 
resolution of a typical neutrino telescope is limited, which leads to an effective averaging of 
any fast oscillation effects. 



4.3 Detector related Systematic Uncertainties 

Solar WIMP searches hold the advantage that background rates can often directly be taken 
from an off-source region, which reduces systematic uncertainties associated with the back- 
ground estimate signiflcantly. Evaluating the signal acceptance is however often a more com- 
plicated undertaking. Since the impact of detector related uncertainties highly dependent on 
the experiment, we can just provide an overview of the most relevant effects. 

Underlying to all detectors is the uncertainty on the neutrino interaction cross section, 
which has uncertainties of a few percent in the GeV - TeV energy region [62]. Somewhat 
larger is the uncertainty for u^{j7^) that is known to about 6%(8%) [35] at the 10 GeV range. 
Neutrino propagation effects through the Earth introduce another source of uncertainty. 

The uncertainty on the energy measurement is another contribution for analyses that 
rely on measurements of the neutrino energy. Considering a heavily instrumented water 
Cherenkov experiments energy resolutions that have been achieved are ±3% for a single 
muon ring and ib(2.5/Y^£'/GeV + 0.5)% for a single electron ring [63]. Further for large 
water /ice Cherenkov detectors also detector specific systematic effects are often dominated 
by photon propagation in the detector medium, limited knowledge of the acceptance and 
absolute efficiency of the individual photomultiplier tubes, as well as reconstruction related 
uncertainties [13, 14]. For vertex contained events, however systematic uncertainties have 
been reduced to about 5% [37]. 
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5 Conclusions 

We have studied the comparison of neutrino fluxes from the Sun with the spin-dependent 
WIMP-nucleon cross section using DarkSUSY. The use of conversion factors for neutrino 
fluxes compared to muon fluxes hold some distinct advantages. Inherently, they do not de- 
pend on muon propagation effects and allow us to treat all neutrino flavors in the same way. 
We have discussed a procedure for the treatment of systematic uncertainties associated with 
the flux conversion and investigated the impact of these uncertainties. While, associated un- 
certainties on the annihilation rate can be signiflcant they are not limiting to measurements. 
Using vertex contained events and spectral event information we have estimated the sen- 
sitivity for neutrino telescope with a corresponding datasets of 200 kton-years and 5 Mton-years. 
For WIMP masses below 100 GeV the sensitivity to the spin-dependent WIMP-nucleon scat- 
tering cross section extends below 10~^ pb. Vertex contained events out perform muon rates 
once the detector size is on the same order as the corresponding track length of the muon 
neutrino induced muon. Even for higher WIMP masses, competitive results can be achieved 
using relatively large opening angles around the Sun that could be used to compensate for 
limited angular resolution. Assuming a comparable angular resolution is achieved the electron 
neutrino channel generally outperforms other flavor channels. 
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